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Low-field magnetotransport measurements on a high mobility (/i=110,000 cm 2 /Vs) two- 
dimensional (2D) electron system on a H-terminated Si(lll) surface reveal a sixfold valley de- 
generacy with a valley splitting <0.1 K. The zero-field resistivity p xx displays strong temperature 
dependence for 0.07< T <25 K as predicted for a system with high degeneracy and large mass. We 
present a method for using the low-field Hall coefficient to probe intervalley momentum transfer 
(valley drag). The relaxation rate is consistent with Fermi liquid theory, but a small residual drag 
as T^O remains unexplained. 
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Two-dimensional electron systems (2DESs) with ad- 
ditional discrete degrees of freedom (e.g. spin, valleys, 
subbands, and multiple charge layers) have attracted re- 
cent interest due to the role of such variables in trans- 
port and in the formation of novel ground states in the 
the quantized Hall regime. In particular, systems with 
conduction band valley degeneracy display a rich param- 
eter space for observing and controlling 2DES behavior 1 . 
Among multi- valley systems, electrons on the (111) sur- 
face of silicon are especially notable because effective 
mass theory predicts the conduction band to be six- 
fold degenerate, yielding a total degeneracy (spin x valley) 
of 12 in the absence of a magnetic field (B). Previ- 
ous investigations of Si(lll) transport using metal-oxide- 
semiconductor field-effect transistors (MOSFETs) with 
peak mobilities p ss 4000 cm 2 /Vs observed a valley de- 
generacy g v of 2 or 6, with the reduced degeneracy at- 
tributed to oxide-induced surface strairP^. 

Here we report transport data on a hydrogen- 
terminated Si(lll) surface (H-Si(lll)) with very high 
mobility (/i = 110, 000 cm 2 /Vs at temperature T=70 mK 
and carrier density n s — 6.7 x 10 11 cm -2 ) with clear six- 
fold valley degeneracy, indicated by the periodicity of 
Shubnikov-de Haas (SdH) oscillations, isotropic low-i? 
transport, and strong T dependence of the longitudinal 
resistivity p xx , consistent with a large gJQ In addition, 
we present a method for using the reduced Hall coeffi- 
cient ru = p X y/(B/en s ) in the B — > limit as a probe 
of valley-valley interactions, using a drag model of in- 
tervalley momentum transfer in multivalley 2DESs. We 
find that the Hall coefficient (and thus, by our model, 
the intervalley drag) becomes strongly suppressed at low 
temperatures (T < 5K); furthermore, although the T de- 
pendence of the drag is roughly quadratic as expected 
from Fermi liquid theory, a small residual drag in the 
T — > limit remains unexplained. 

To create and probe a high- mobility electron system on 
a bare surface, we fabricate a device similar to a four wire 
MOSFET, with the critical difference that we replace the 
Si-SiC>2 interface with a H-Si(lll) surface adjacent to a 



vacuum cavity^ 1 . The main processing enhancements in 
the device discussed here relative to our prior samples 
are a higher resistivity Si(lll) substrate (p ~ 10 kf2-cm) 
and final H-termination and bonding performed in an 
oxygen-free (< 1 ppm) environmental. The resulting de- 
vice has a very high mobility which, as Fig. [I] illustrates, 
increases monotonically with n s and thus is likely limited 
by charged impurity scattering due to residual surface 
charge. The sample was probed via standard Van der 
Pauw measurements (Fig. [2] inset) using low frequency 
(7-23 Hz) lock-in techniques in both He-3 and dilution 
refrigerators. 

At B=0, p xx is strongly affected by T (Fig. [ljj) and 
at low T displays a metal-insulator crossover near n M1T — 
0.9 x 10 11 cm -2 . Above this density, the device is clearly 
metallic, with p decreasing by a factor of 3-4 from 5 
K to 70 mK at n s = 7 x 10 n cm~ 2 . Compared with 
Si (100), Si(lll) has a larger density of states effective 
mass {m*ii — 0.358m e vs. ml m = 0.190m e ), and a 
larger g v (6 vs 2) would lead to a much larger density 
of states at the Fermi level. Consequently, electrons 
on Si (111) should display much stronger screening and 
therefore a stronger T dependence of p x J^. Our observa- 
tions appear to be qualitatively consistent with such pre- 
dictions. At intermediate densities, p xx is non-monotonic 
in T, similar to behavior discussed elsewhere^. The data 
presented in the remainder of this paper was taken at a 
fixed density of n s = 6.7 x 10 11 cm~ 2 . 

The low-B SdH oscillations of the present sample re- 
veal minima at filling factors v — hn s /eB spaced by 
Av = 12 below 0.5 T and Av = 6 above this point (Fig. 
[2J. Note, however, that the 12-fold-degenerate minima in 
p xx occur at odd multiples of 6 (54, 66, 78, 90...) rather 
than even as would be expected if the effective mass m* 
and g-factor g* are equal to their band values. This may 
indicate a independent valley splitting and/or enhanced 
Zeeman splitting at low fields (though unlike the spin- 
dominated gaps seen in Si(100^and Si(lll)0MOSFETs, 
our observations persist far into the metallic regime). At 
B > 2 T valley degeneracy lifts, eventually resulting in 
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FIG. 1: a) Hall mobility vs. carrier density n s . Up- 
per inset: Sixfold valley structure of Si(lll). Lower inset: 
Schematic device cross-section showing two of four contacts. 
The (500/im) 2 vacuum cavity created by contact bonding pre- 
serves the H-Si(lll) surface and serves as the gate dielec- 
tric, b) p vs. T for n s ranging from 0.65 (■ at top) to 
7.03 x 10 11 cm -2 (♦ on bottom). The crossover from 'metallic' 
(dp/dT > 0) to 'insulating' (dp/dT < 0) behavior occurs near 



n MIT = 0.9 x 10 1 



integer quantum Hall features appearing at intervals of 
Av = 1 from v = 10 to 1. 

Another effect of the weak valley splitting is the obser- 
vation of isotropic transport (R xx — R yy to within « 7%) 
as shown in Fig. [2] (inset). If charge is evenly distributed 
among a suitably symmetric set of valleys (discussed be- 
low), the total resistivity becomes isotropic. A lifting 
of the valley degeneracy can shift the valley population 
into an asymmetric distribution, causing anisotropy in p 
and other transport effects^. While previous H-Si(lll) 
devices have shown this anisotropy^, the absence of such 
asymmetry in the device presented here suggests that the 
valley splitting is quite small at B=0. When a small B is 
applied, the baseline R xx and R yy values remain similar 
but the SdH response is different: the p yy minima appear 
out of phase with the p xx minima while retaining the es- 
sentially 12- fold periodicity (Fig. |2j. This may result 
from a slight population imbalance due to a very small 
valley splitting, which may increase further for B > 0. 

In order to provide a quantitative bound on the zero- 
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FIG. 2: Low-B SdH oscillations for orthogonal current di- 
rections. Below 0.5 T minima appear every 12 levels, indicat- 
ing a twofold spin and sixfold valley degeneracy. Oscillations 
continue as high as v w 234 (B « 0.1T). Above 0.5 T fur- 
ther structure becomes evident and minima appear in steps 
of Au — 6. While this interval size is essentially the same 
for both traces, the phase is clearly different. Inset: Base T 
anisotropy ratio R xx /R yy at B — 0; diagram shows current 
configurations for 4- wire measurements of R xx , Ryy, and R xy . 

field valley splitting, we measure the intrinsic (i.e. B- 
indepcndcnt) level broadening, which is normally charac- 
terized in terms of the Dingle temperature Tjj, or equiva- 
lently by the quantum lifetime r q = h/2nkBTr) 1 by exam- 
ining the evolution of the SdH oscillations as a function 
of T. If the energy level spacing is E gap = SB (for some 
constant S) the amplitude of these oscillations should 
be R e- 2n2kBTD / SB t/smh(£) where £ = 2-K 2 kT/SBW^ 
Thus from the T dependence of the amplitudes (unaf- 
fected by the phase anomalies noted above) we can de- 
termine T q as well as the B-dependent gap size 5. 

From the p xx oscillations we obtain a value of 8 = 
(2.69±0.11) K/T (Fig. [3). Because our gaps occur at odd 
multiples of six, the simplest analysis would treat them 
as spin gaps, which would correspond to an enhanced 
g* = 4.0 ± 0.2. However, the anisotropy that emerges 
when B > suggests B-dependent valley splitting is also 
present, confounding a simple interpretation of 5. Re- 
gardless of the gaps' origin, from 5 we can compute the 
quantum lifetime r q ss 12 ps, which is quite close to the 
transport lifetime tq ~ 18 ps obtained from p xx (B = 0) 
as described below; this corresponds to a Trj (and thus 
an upper bound on the valley splitting) of 0.1 K. 

Having established the sixfold valley degeneracy of our 
2DES at B=0, we now turn to the role of this degener- 
acy in carrier scattering; in particular we consider the 
effect of momentum exchange between valleys on low B 
transport. Semiclassical transport in multiple anisotropic 
valleys can result in additional (non-oscillatory) B depen- 
dence in both p xx and p xy at low fields (seen, for example, 
in the overall positive slope of the data in Fig. |3]P. This 
low B behavior can provide information about valley- 
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FIG. 3: Low-B SdH oscillations at several temperatures. 
We compute the size AR of the oscillations via linear in- 
terpolation of the minima and maxima across a wide range 
of t^. We then plot AR/T vs T to extract the gap size 
5 = E gap /B = 2.69 K/T. Inset: Using this S, we plot 
AR/(£/ sinh(f)) vs. 1/B to compute the Dingle temperature 
Td = 0.1 K. Note how the data collapse to a single line. 

valley interaction effects. To see this, we first consider 
the case of non-interacting valleys that are identical up 
to rotations Z(9) in the x-y plane for some 9 < it that 
defines the rotational symmetry of the whole set. The 
Drude resistivity for a single valley (with proportional 
density n s /g v ) for coordinates aligned to the symmetry 
axes of the valley is given by 
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where To is the transport lifetime associated with momen- 
tum transfers from the 2DES to the lattice (both intra- 
and inter- valley scattering). The resistivity of the j th 
valley is pj = Z(j9)p Z(—j9) and the total p will be 
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where m = (mi + m2)/2, uj c = eB/m* is the cyclotron 
frequency, m* — y/m\m2, and we define $ = (m/m*) 2 . 
For u> c tq 1 both p xx and p xy are given by their respec- 
tive classical values p xx = fh/e 2 n s TQ and p xy = B/en s . 
At B — however, both are suppressed by the factor 

< 1, with equality only in the case of isotropic val- 
leys (mi = m-i). In the case of Si(lll) (using the band 
masses mi = 0.190m e , m 2 = 0.674m e , m* = 0.358mJ^) 
we find i = 0.686. 

The preceding discussion treats the valleys as indepen- 
dent channels. Strong valley- valley coupling, as might 
arise from Coulomb interactions between electrons, will 
tend to suppress this correction as all electrons move in 



concert. To make this idea more rigorous, we model in- 
tervalley effects as a drag interaction between valleys that 
conserves total 2DES momentum while damping the rel- 
ative momenta between valleys. This is distinct from in- 
tcrvalley scattering probed via weak localization^ 3 which 
requires a short-range interaction potential that does not 
conserve 2DE S mome ntum. Following the kinetic ap- 
proach used ir i 14 ! 15 * 16 ! for multi-band systems we obtain 
a set of coupled equations: 
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Here Mj is the mass tensor of the j valley. 



Mr, 1 = 



MJ 1 + M, 
and 



, is the reduced mass tensor of the j-k system, 
t vv is the drag relaxation time (assumed constant 
and isotropic). Combining opposite valleys (which have 
the same Mj), we have three valley pairs. Substituting 
jfc = ^ifcevk = n s eVk/3 we then solve the equation E = 
pJ = Pih+h+h) for p. This gives 
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where A = 6det(M 7 - fe )/det(M J ) = 6/(3$ + 1). In the 

absence of intervalley interaction, ► and we recover 

Eq. Conversely, when t vv <C tq we effectively wash 
out the multi-valley correction. Thus by measuring p xx 
and p xy in the B = limit we can solve for tJ 1 and t^, 1 . 

Figure E] (left axis) shows such a measurement of 
Tji = p xy J[B / en s ) vs. T, averaging results from orthog- 
onal directions (Fig. |2| to remove mixing from p xx and 
Pyy. We determine the density n s = 6.7 x 10 11 cm -2 
from R xx minima at v — 18 and v — 6 and find it to be 
insensitive to T. B was held fixed at ±50 mT while T 
was swept both up and down to ensure consistency. Tak- 
ing a slope from these points gives a measure of rjj near 
B~0. Above 5 K, th is very close to its classical value, 
while below 5Krg drops rapidly with T before settling 
to a value of 0.65 at T = 90 mK. Interestingly, 0.65 is less 
than the lower bound of 0.686 predicted for the drag-free 
limit. Because the measurement is based on data taken 
at B ^ 0, we expect this to be an overestimate of r#, 
especially at low T where the uj c tq terms in Eqs. Q and 
([5]) are largest. The simplest adjustment we can make to 
our model to incorporate this discrepancy is to allow 
to approach a constant negative value at low T. 

On the right axis of Fig. |4]we plot the T dependence of 
the extracted t~ v 1 +(364: ps) -1 (offsetting t~J by the base 
temperature value to remove the divergence) as well as 
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FIG. 4: Left axis: Reduced Hall coefficient rH measured 
near B — {—)■ th = 0.686 is the lower bound predicted 
by the simplest model. Right axis: Lattice scattering rate 

Tq -1 (. . .) and valley drag relaxation rate t~J ( ) vs. T, 

computed from th and p xx data via Eqs. Q and The 
inverse quantum lifetime r^ 1 (-fa) is calculated from low-T 

SdH oscillations. The shows the damping rate expected 

for a Fermi liquid (see text). 

the lattice scattering rate Tq 1 . The dashed magenta line 
plots the T-dependent electron-electron (e-e) scattering 
rate theoretically expected^ for a single- valley 2D Fermi 
liquid t~1 ~ "if (^:) 2 i with a pref actor of 1.9 determined 
by fitting. Although earlier wor k has id entified the sen- 
sitivity of Th to e-e interaction ;] 14 * 18 * 1 -^ such corrections 
are quite small (~l-4%) at high densities (n s ^> n MIT ). 
Furthermore, these models predict rjj — ► 1 in the T — > 
limit, whereas our model treats the multi-valley effects 
of Eq. ^ as intrinsic at T = 0. 

Several open questions remain regarding this data. 



First is the dominance of the odd gaps and the low-field 
phase anisotropy in the SdH data. Although the exper- 
imental setup did not allow for tilted B field measure- 
ments, such experiments could help identify the roles of 
cyclotron, spin, and valley splitting in forming these gaps. 
Second, we note the disparity between the measured th 
at base temperature and the lower bound given by our 
model, which we have presented in terms of negative drag 
for mathematical simplicity. Known corrections to p xv , 
for instance due to disorder, would increase ^^EISEI 
not decrease it as reported here. Negative drag has been 
discussed as a possibl e conse quence of electron correla- 
tion in bilayer systems ^ * 21 * 22 !. Small sample anisotropics 
(neglected in our model) could also play a role. Another 
possibility is that anisotropic enhancement of m* could 
change the T = limit of Eq. Q (note that isotropic 
enhancement would not change rjy, which depends only 
on mass ratios). Alternatively, a more sophisticated the- 
ory of valley- valley interactions could modify the T — * 
limit of our simple model. 

Finally, we consider why sixfold degeneracy persists 
in this device given the ease with which mechanisms 
such as mi sorient ation, disorder, and strain can lift this 
degenerac y^ * 10 * 23 ! and that previous work on H-Si(lll) 
found the g v = 6 ground state split into a low energy 
g v — 2 band with a g v = 4 band ~7K above itP. If that 
gap was in fact produced or enhanced by surface disor- 
der, perhaps the higher mobility of the present sample 
can account for the difference. Further work on the re- 
lationships between surface preparation, device mobility, 
and valley splitting is presently underway. 
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